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SUMMARY 

Tile rotating— stall ciiaracteristics of a 0,9 Lnb— tip ratio rotor were 
investigated. Stall patterns consisting of two^ th.ree^ and one total- 
span stall zones developed in that order upon reduction of flow coeffi- 
cient. The one -stall -zone pattern caused the most severe pressure, tem- 
perature, and flow fluctmtions . 

As the flow decreased within the stall zone, detailed measurements 
showed an increase in static pressure within the stall zone upstream of 
the rotor and a static-pressiore decrease downstream of the rotor. These 
pressure fluctiiations were of sufficient amplitude in the case of the 
one- stall- zone pattern to result in a pressure drop across the rotor 
within the stall zone. 

All stall zones rotated in the direction of rotor rotation relative 
to the compressor casing and at a speed proportional to hut less than 
rotor speed. The measured stall- propagation rates were compared with 
those predicted hy theory. The relative stall-propagation rate increased 
with increasing stall- zone size. 


INTRODUCTION 

Foiir extensive experimental investigations reported in references 1 
to 4 show that stall in axial-flow and centrifugal compressors may mani- 
fest itself in asymmetric flows. This type of stalled condition has come 
to be known as propagating or rotating stall (ref. 4 refers to the phe- 
nomenon as "type 1 surge"). The asymmetric flows consist of one or more 
low-flow regions encompassing several blade passages. These low- flow 
regions rotate about the compressor at a speed proportional to but less 
than rotor speed and in the direction of rotor rotation relative to an 
absolute frame of reference. 
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When rotating stall is encountered, an abrupt decrease in compressor 
pressure rise may occur when a single total-span stall zone develops. 

This type of stall is often called abrupt stall. The compressor then 
operates on a stalled branch of the compressor performance characteristic 
as shown in figure l(a). Upon increasing the weight flow, the compressor 
recovers from the severe stalled condition at a higher weight flow than 
that where the drop in pressure originally occurred, thus forming the 
hysteresis loop indicated in figure l(a). Reference 4 reports two stalled 
branches for one compressor^ however, the data indicate that a compressor 
usually has one preferred stalled operating characteristic at any given 
speed. References 1 to 4 report some compressor configurations where a 
continuous performance characteristic similar to that shown in figure 
1(b) is obtained when rotating stall is present. Generally? in these 
cases only a portion of the blade span is stalled. This type of stall 
is often called partial-span or progressive stall. 

In addition to the effect of rotating stall upon performance, ref- 
erence 5 suggests that rotating stall is a source of resonant excitation 
that may restilt in blade failure. Compressor blade failures directly 
attributable to rotating stall are reported in reference 6. 

The flow fluctuations associated with rotating stall are distinctly 
different from the phenomenon generally called "surge." During operation 
with rotating stall, the net flow rate throtigh the compressor and the 
torque necessary to drive the compressor remain essentially constant, 
whereas these quantities fluctuate during surge. 

The problem of rotating stall in a single blade row is analyzed in 
references 7 to 10. In references 7 to 9 the blade row is considered to 
have an infinitesimal chord length (actuator sheet), whereas the effect 
of blade chord is included in reference 10. The analyses show qijalita- 
tively that flow patterns of the rotating- stall type may occTir in axi- 
symmetric rigid blade rows with steady uniform inflow to the blade row. 

To obtain a better understanding of the stall phenomenon, an in- 
vestigation of the rotating- stall characteristics of a 0.9 hub- tip ratio 
rotor has been carried out at the MCA Lewis laboratory. The high hub- 
tip ratio rotor is typical of an exit stage of a multistage compressor 
and corresponds approximately with the two-dimensional single-blade-row 
models considered in the theoretical treatments. The investigation was 
made in order to provide more detailed information regarding the stall 
mechanism, pressure, temperature, and flow fluctuations and their axial 
extent, and the phasing relation of the fluctuations. Results are com- 
pared with theoretical analyses where applicable. 
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APPARATUS AND INSTRUMENTATION 
Rotor Design and Installation 

The 14 -inch-diameter rotor used in this investigation has 31 alumi- 
num blades of 0.9 huh-tip ratio. The blades are untwisted and untapered. 
NACA 65(12) -10 blower blade profiles were used. The rotor blade setting 
angle was 28.5^^ and the inlet-air angle for minimum losses as determined 
from cascade tests was 45^. Figure 2 shows the design vector diagram. 

The rotor test setup is shown schematically in figure 3. Air en- 
tered through the orifice tank and a motor-driven inlet throttle into 
an inlet tank where screens were used to straighten the air at the com- 
pressor inlet. The compressor air discharged into a receiver and then 
through an outlet throttle to the laboratory altitude exhaust system. 


Flow-Fluctuation Instrumentation 

A const ant -teii 5 )erature hot -wire -anemometer system such as reported 
in references 11 and 12 was used to measure flow fluctuations. The wire 
element of the anemometer probe was made of 0.0002 -inch-diameter tungsten 
wire mounted parallel to the probe ax:is and oriented in the stream so 
that the wire element was not in the wake of its supports. The anemometer 
output signal was filtered electronically to pass up to 2000 cps and re- 
corded photographically from a dual -beam direct -coupled cathode -ray 
oscilloscope. 

Hot -wire -anemometer probes were insta3JLed at axial stations 2, 4, 

5^ 6^ 7, and 8 (fig. 4). Probes were installed at several angular 
spacing s at station 4. 


Pressure -Fluctuation Instrument at ion 

Pressure fluctuations during stall were measured with a calibrated 
inductance -type pressure pickup of the kind reported in reference 13. 

For the purpose of the investigation^ the pressure pickup was installed 
in the compressor casing as shown in figure 5. The pressure pickups were 
inserted into the compressor casing with the static tap flush with the 
inner surface of the con 5 )ressor casing at axial stations 3 and 7 
(fig. 4). Radial measurements were obtained from "extension disks" from 
the static tap. The amplitude response of the pickup was "flat" to a 
frequency of 1000 cps (ref. 13), which was well above the stall frequencies 
encountered. The output of the pickup was photographically recorded from 
a dual -beam direct -coupled cathode-ray oscilloscope. 
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Temperature -Measurement Instrtimentation 

The average temperature rise across the rotor row was measured by 
Iron-constantan thermocouples mounted in stagnation shields at axial sta- 
tions 4 and 7. The temperature fluctuations during rotating stall were 
measured at axial stations 4 and 7 by a compensated thermocouple made of 
chromel-constantan 0.003 inch in diameter and also with a 0.0002- inch- 
diameter tungsten-wire resistance thermometer. Both of these devices are 
discussed in reference 14. The re si stance -thermometer and compensated- 
thermocouple outputs were viewed on a dual -beam oscilloscope. 


PROCEDURE 

Rotor Performance Characteristic 

The steady-state performance of the rotor is expressed in terms of 
pressure and flow coefficients. The curve showing the relation between 
these two parameters is referred to as the static performance character- 
istic. The pressure and flow coefficients, as defined in the symbol list 
of appendix A, were determined from data from 35 static -pres sure taps 
along the inner and outer casings and from the weight flow and tempera- 
ture. The flow coefficient was varied from choked flow to a low-flow 
condition by means of the outlet throttle. 

Investigation over several speeds showed that the stall zones rotated 
at a speed proportional to the rotor speed. The data reported herein were 
taken at = 8000 rpm (rotor tip speed of 500 ft/sec). An absolute 

inlet tank pressure of 25 inches of mercury was maintained at all weight 
flows . 


Measurement of Flow, Pressure, and Temperature Fluctuations 

The voltage output of the hot-wire anemometer was converted into 
mass -flow fluctuations by the method of appendix B. The number of stall 
zones was found by displacing hot-wire probes around the circumference 
at several different angular spacings at axial station 4 and applying the 
analysis of appendix B. The stall frequency was determined by forming 
Lissajous figures on the oscilloscope with an audio-frequency oscillator. 

The amplitudes of the static-pressure and total-temperature fluctua- 
tions were measured at axial stations 4 and 7 (fig. 4) for each stall 
pattern obtained, along with the flow fluctuations at axial stations 2, 

4, 5, 6, 7, and 8. 
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RESULTS AND DISCUSSION 
Performance and Stall Characteristics 

The static performance characteristics for the annulus between sta- 
tions 1 and 5, for the rotor between stations 5 and 6, and for the com- 
plete unit between stations 1 and 6 are given in terms of pressixre and 
flow coefficient in figure 6. While the pressure coefficients for the 
inlet annulus and for the complete unit continue to rise upon reduction 
of flow coefficient, the pressure coefficient for the rotor begins to 
fall off at approximately <P = 0.9, until at cp = 0.75 a rotating- stall 
pattern of two zones develops. At this point the over-all performance 
as well as the rotor performance is affected by the rotating- stall condi- 
tion. However, the amplitude of the flow fluctuation is not of suffi- 
cient magnitude to affect the pressure coefficient of the inlet ann\olus. 
Upon further reduction in flow coefficient, a three-stall-zone pattern 
develops at (p = 0.69 and persists to a flow coefficient of 0.64, where 
a single stall zone covering about one -third of the annulus and producing 
large flow fluctuations develops. The single-zone stall pattern has a 
considerable effect on the inlet -annulus and over-all pressure coeffi- 
cients. The single rotating- stall zone remained throiigh a flow coeffi- 
cient of 0.52. The valve arrangement of the system did not permit flow 
coefficients between zero and 0.52 to be obtained conveniently, so no 
data were teiken in this flow-coefficient range . 

The rotor was run for a short time with the valve completely closed 
to determine the rotating- stall characteristics at zero weight flow. 

One stall zone was found rotating at half rotor speed. 

The rotating- stall patterns are summarized as follows for a rota- 
tional speed N of 136 rps: 


Range of flow 
coefficient, 

<P 

Nvnnber of 
stall zones, 
X 

Frequency of 
stall zone 
passing hot- 
wire probe, 

^s> 

cps 

Rotational 
speed of 
stall zones, 

rps 

Ratio of 
stall -zone 
speed to 
rotor 
speedy 
h/N 

0.75 - 

0.69 

2 

175 

87.5 

0.64 

.69 - 

.64 

3 

21 e, 

92.0 

.68 

.64 - 

.52 

1 

60 

60.0 

.44 

0 


1 

68 

68.0 

.50 


If the rotor were operated at the flow coefficient of 0.52 for 
several minutes, the stall pattern would change from X = 1 to X= 4. 
The four stall zones mi^t then later regroup into the original one zone. 


6 


NACA TN 3518 


This breaking down into four zones would recur on subsequent operation 
at this point, indicating that the one stall zone was not unique at this 
point even though it was the first to appear. 


Rotor Pressure-Rise Coefficient 


The two-dimensional pressure characteristic of the rotor at midspan 
(r = 0.95) in terms of cascade parameters is given in figure 7. The 
pressure-rise coefficient ^v/%^ increases until an inlet angle of 50° 

is reached and then remains constant until approximately 52°, where ro- 
tating stall is encountered and the pressure-rise coefficient decreases. 
Because the flow is fluctuating with time, the pressure rise measured 
during operation with rotating stall is a time -averaged value and not 
necessarily the true static characteristic. Since rotating stall pro- 
duces alternating high- and low-flow areas, the inlet angle also oscil- 
lates over a wide range. 


The rotor pressxxre rise obtainable assuming two-dimensional, isen- 

1 ^ 

tropic, compressible flow and a constant exit angle of 221 is also given 


in figure 7. The slope of the curves is the same except \diere the inlet 
angle deviates more than 5° from the design angle. The rotor pressure - 
rise coefficient is approximately 75 percent of the isentropic value at 
the design inlet angle. The rotor blade choked at an inlet angle of ap- 
proximately 36° and a relative inlet Mach number of about 0.73, as is 
indicated by the sharp reduction in pressure-rise coefficient. 


o 


Flow Fluctuations due to Rotating Stall 

Variation in amplitude of flow fluctuations with axial position . - 
Flow fluctuations caused by rotating stall are shown in figure 8. The 
average flow fluctuations are reported at the mean radius (r = 0.95). 

In general, the anemometer signal obtained at station 6 was too "hashy" 
to analyze quantitatively for flow- fluctuation amplitudes, and no data 
are reported for this station. As may be seen from figure 8, the flow 
fluctuations for X = 1 were the most severe and of nearly the same 
amplitude at n.l 1 axial measuring stations. The fluctuations for two- 
and three-zone stall patterns were largest near the rotor and diminished 
in amplitude both upstream and downstream of the rotor. 

Representative oscillograms from axial surveys are reproduced in 
figure 9. Multiple traces appear on some of the oscillograms, because 
the oscilloscope beam swept more than once during the film exposure time. 
All the traces appear in pairs with one anemometer trace from station 4, 
selected for comparison. The oscillograms in figures 9(a) and (b) for 
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the two- and three -stall patterns^ respectively^ show clearly the small- 
amplitude flow fluctuations upstream and the relatively larger fluctua- 
tions in the downstream annulus. Where one stall existed (fig. 9(c)), 
the flow fluctuations were of nearly the same amplitude at all stations, 
with little change in wave shape downstream of station 2. 

Variation in amplitude of flow fluctuation along hlade span . - As 
might be expected for a short blade height, flow fluctuations caused by 
rotating stall existed over the entire blade span for all rotating- stall 
patterns. The results of a radial suivey at station 5 immediately pre- 
ceding the rotor are given in flginre 10. The three-stall-zone pattern 
exhibited the greatest radial variation in fluctuation amplitude , The 
fluctuations become more equally distributed spanwise as the mmiber of 
stall zones is reduced. The flow fluctuations with one stall zone are 
uniformly distributed over the entire span and as pointed out earlier 
are of much greater amplitude than either the two- or three -stall-zone 
flow fluctuations. 

Typical oscillograms from the radial surveys at station 5 are given 
in figure 11 for the two-, three-, and one-stall patterns. There is a 
pronounced change in wave shape between the two- and three-stall-zone 
patterns of figures 11(a) and (b) and that for the one -stall pattern in 
figure 11(c) . 

It seems altogether likely that some reversed flow took place near 
the rotor with the single-zone stall pattern. The upswing within the 
stall-zone "blip" on figure 11(c) may be an indication of reversed flow. 
Reversed flow or back flow within the stall zone is sometimes observed 
(e*g*^ ref. 4). Such observations are usually made by detecting the 
motion of tufts mounted upstream or downstream of the compressor with 
stroboscopic light. Since the anemometer probes used in this investiga- 
tion are insensitive to flow direction, reversed flow could not be posi- 
tively identified. 


Pressxrre Fluctuations due to Rotating Stall 

Radial surveys with the pressure pickup at stations 3 and 7 indicated 
that the pressure variations during stall were essentially constant over 
the radius for all stall patterns. Typical oscillograms of the pressure 
fluctuations are shown in figure 12. For all stall patterns, the ampli- 
tude of the pressiire fluctuations downstream of the rotor was approximate- 
ly 60 percent of the upstream fluctuations. The hipest -amplitude fluctua- 
tions occurred for the single-zone stall pattern (fig. 12(c)). The 
pressure fluctuations at stations 3 and 7 were nearly 180° out of phase, 
as shown by the oscillogram in figure 13, which was obtained during opera- 
tion with the single -zone stall pattern. 
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The pressiire fluctuations at stations 3 and 7 are shown schemati- 
cally in figure 14. The two- and three -zone stall patterns produced ap- 
proximately a ±10-percent variation in pressure rise between stations 3 
and 7. The single -zone stall pattern^ however, resulted in a much 
greater variation, as shown in figure 14(c). The maximum static pres- 
sure at station 3 d-urlng operation with the single-zone stall pattern 
was sli^tly greater than the inlet total pressure (station l) j the minl- 
Ttrum press\ire at station 7 was only slightly greater than the average pres- 
sure at station 3. That is, within the stall zone the static pressure at 
station 3 was increased to approximately the inlet stagnation pressure, 
while the static pressure at station 7 decreased to a value slightly 
greater than the average static pressure at station 3. Within the stall 
zone there was a pressure drop across the rotor. This drop is shown in 
terms of pressure coefficient in figure 15, where the instantaneous pres- 
sixre coefficient between stations 3 and 7 is plotted against time. 

The hot -wire -anemometer data obtained during operation with the 
single-zone rotatlng-stall pattern indicated that the flow coefficient 
at station 4 varied from nearly zero in the stall zone to a value between 
0.85 and 0.92 between stall zones. The uncertainties Involved in inter- 
preting the anemometer data prevent a more accurate determination of the 
flow coefficient between stall zones. The steady- state or static per- 
formance of the rotor and inlet annulus is shown in figure 16 in terms 
of pressure coefficient for comparison with the instantaneous pressure 
coefficients obtained with the single-zone stall pattern. According to 
figure 14(c), the maximum pressure coefficient across the rotor '^3_y 

was 0.805 and the minimum press\rre coefficient in the inlet annulus ^ 

was -0.784. These values correspond to a flow coefficient <P of 0.87 on 
the static performance ctirves in figure 16, compared with a flow coeffi- 
cient between 0.85 and 0.92 estimated from the anemometer data. The maxi- 
mum pressure coefficient outside the stall zone for the single-zone 

stall pattern is therefore approximately the same as the maximum steady- 
state (unstalled) pressvire coefficient. 

Inasmuch as the flow fluctuations obtained with the two- and three - 
zone stall patterns were somewhat smaller than those with the single-zone 
stall pattern, the fluctuations in pressure coefficient were correspond- 
ingly less, as indicated in figures 14(a) and (b) . 

The data indicate that the pressixre coefficients between stall zones 
correspond approximately to the steady- state values at the steady- state 
flow rate. A similar comparison for conditions within the stall zone 
cannot be made, because steady- state data at low (stalled) flows were not 
obtainable . 
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Temperature Fluctuations due to Rotating Stall 


The measured average temperature rise across the rotor is plotted 
against flow coefficient in figure 17. The temperature rise increases 
rapidly at low flow coefficients where rotating stall is encountered. 

The amplitude of the temperature fluctuations at station 7 during rotat- 
ing stall was 4° F for the two- and three-stall patterns and 25° F for 
the single-zone stall pattern. 

Figure 18 represents the approximate fluctuation in total tempera- 
ture at station 7 diiring operation at a flow coefficient of 0.63 with 
the single -zone stall pattern. The average temperature rise under these 
conditions is 23® F (fig. 17) . The minimum temperature rise occurring 
in the flow area outside the stall zone is 16° F^ and the maximum rise 
in the stall zone is 41° F. 

The analysis of the pressure oscillograms and anemometer signals 
indicated the through-flow velocity in the stall zone to be essentially 
zero. Assuming^ then, that the air in the blade passages lying within 
the stall zone rotates at blade speed, the temperature rise across the 
rotor in the stall zone would be 


which is in good agreement with the measured value of 41° F. The quasi - 
steady-state analysis indicated a flow coefficient of 0.87 outside the 
stall zone. Under steady flow conditions at this flow coefficient, a 
temperature rise of 18° F was measured (fig. 17), which is close to the 
16 F temperature rise measured at the same flow coefficient under un- 
steady flow conditions. 


Discussion of Rotating Stall and Comparison with Theory 

The mechanism by which the stall zones propagate from blade to blade 
is usually explained in a qualitative manner as follows. Consider an 
infinite row of blades operating at the stall point in two-dimensional 
flow (see fig. 19). Now suppose that a small pertijrbation in the incom- 
ing flow causes blade 2 to stall and consequently become incapable of 
supporting the existing pressure rise across the blade row. As a con- 
sequence of stall, the flow will be reduced through the separated blade 
passage. The incoming flow will be diverted around the local low- flow 
area, and hence the angle of attack on blade 3 will increase while the 
angle of attack on blade 1 will decrease. Blade 3 will stall, and the 
low-flow or stall zone will propagate in the direction of the component 
of velocity parallel to the cascade. After a fraction of a second, the 


H 


O 



2 


= 38° F 
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stall-zone propagation rate and size would presimably reach some equilib- 
rium value covering several blade passages. For an annular cascade^ 
the stall zone would propagate around the aixis of the cascade and is thus 
called rotating stall. 

The rotor used in this investigation corresponds approximately to 
this model. The hub-tip ratio is 0.9, so that the flow is nearly two- 
dimensional. Relative to the rotor, the stall zones propagate in the 
direction of the tangential velocity but at a velocity less than the 
rotor speed, so that in an absolute sense the stall zones rotate in the 
direction of the rotor rotation. 

The qualitative explanation of stall propagation gives the correct 
direction of rotation of the stall zones and also implies the increase 
in static pressure within the stall zone ahead of the blade row. This 
explanation also implies that operation at the stall point is unstable. 

If the flow were stable, the perturbation would be damped and steady 
rotating stall woijld not be obtained. 

A prediction or estimation of the rotative speed of the stall pat- 
tern and the number of zones formed requires additional considerations. 
References 7 to 10 present linearized analyses that result in expressions 
for the stall-propagation rate. There are no methods available for esti- 
mating the niimber of stall zones that will form. The theories consider 
two different forms of a single blade row. Marble (ref. 9) and Sears 
(refs. 7 and 8) represent the blade row by an infinite number of small- 
chord blades, so that the configuration they use is in effect an actuator 
sheet. Stennlng (ref. 10) considers a cascade of blades of finite chord 
length . 

Sears considers a flow represented by an Irrotational flow field 
entering the rotor, an actuator sheet acting as the rotor, and a rota- 
tional field leaving the rotor. The vorticity following the actuator 
sheet is determined either from consideration of the shedding of vortices 
by the blades as they pass through the stall zones ("airfoil theory") or 
from the variation in downstream total pressure ("channel theory") . 

In the airfoil theory, expressions are derived for the blade circula- 
tion in terms of the velocity distribution. The induced inflow velocities 
at the actuator sheet are calculated from the shed vorticity. After in- 
troducing a phase angle A between lift and angle of attack, the condi- 
tions required for the blade -loading variations to be just supported by 
the induced velocities are determined. The conditions required are as 
follows ; 
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' ^ ( 1 tan^ P - 1 \ //-, tan^ P - 1 2 

N ~ \ ^ Z tan p tan A/ y\ 2 tan p tan A/ ^ tan p tan A (2) 

In order for stall propagation to occur with the slope of the lift c\irve 
m just slightly negative, the phase shift A must approach 90°. The 
experimental stall -propagation rates are compared with those obtained 
considering the negative root of equation (2) in figure 20. The phase - 
shift values A required for close agreement with the data are 80° and 
^ 85°, as seen from figijre 20. At zero weight flow (p = 90°), the theory 

g would predict either no stall or a stall propagating at rotor speed, 

O) since h/N is zero. Experimentally, a value of h/N = 0.50 was ob- 

tained. 

The channel theory (ref. 8) utilizes relations derived from cas- 
cade tests - namely, that the relative outlet-air angle and the pressure 
rise across the blade row are functions of the inlet-air angle. The in- 
duced velocities are calculated from the variation in total pressure 
downstream of the rotor. A phase angle 6, similar to that used in the 
airfoil theory, is introduced; and the procedure used in the analysis is 
much the same as for the airfoil theory. The conditions required for 
the downstream total-pressure variation to be just supported by the in- 
duced velocities are as follows: 


/ /h\2 


+ \n) (sec^p) 

( 3 ) 

= -cot p cot (6 + 2P) 

( 4 ) 


A plot of equation (4) for 6 of 0°, 10°, and 20° is shown in figure 
21. (The 6 = 0° curve will be discussed f\orther in connection with 
Marble's theory.) Values of 6 of 10° and 20° bring the channel theory 
in reasonable agreement with all the data except at a flow coefficient 
of zero, in which case a phase shift of zero is required for agreement 
between theory and experiment. 

Marble in reference 9 like Sears in reference 8 represents the blade 
row as an actuator sheet. The press\jre-rise curve across the disk is as- 
sumed to be a linear function of the inlet-air angle for unstalled flow 
and to drop discontinuously to zero for stalled flow. Pressure is used 
as the dependent variable in the analysis instead of the velocity com- 
ponents used by Sears. The asymmetric flow is introduced by considering 
a stall zone of variable circumferential extent to exist in the actuator 
disk. Two cases are considered: a "simple cascade" solution, where the 

deviation angle remains constant for all values of incidence angle, a n d 
a "general cascade characteristics" solution, where the cascade pressure 
rise and the deviation angle are considered functions of the incidence 
angle. A phase angle between angle of attack and pressxnre rise is not 
considered. If the phase angle 6 is assumed zero in Sears' channel 
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theory^ the expressions for stall-propagation rate obtained by Marble for 
the simple cascade and by Sears are identical. This relation is shown 
in figure 21. In addition to propagation rate. Marble finds an expres- 
sion for stall-zone width. 

The stall analysis by Stenning (ref. 10) considers a cascade of 
finite chord length. It is assumed that the effective outlet area of the 
blade passage is reduced (flow coefficient reduced) as a result of stall 
and that the flow coefficient of the blade passages is determined by the 
inlet angle. The static pressure is assumed constant downstream of the 
cascade in spite of the asymmetric flow pattern. The data obtained in 
the present investigation indicate that the downstream static pressure is 
not a constant. The effect of the variation in downstream static pressure 
on propagation rate has, however, not been ascertained. Stenning obtains 
a solution for the stall-propagation rate from momentum considerations 
within the cascade. The result obtained where the blade row is assumed 
to consist of a large number of short -chord blades is 

, h _ cot^p 
^ cos^y 

One of the results of Stenning 's work considering the finite- chord- 
length blade row is that the relative stall-propagation rate increases 
with the size of the stall zone . The experimental results of the present 
investigation agree with this conclusion, as shown in flgvn*e 22. In this 
figiire the relative propagation rate is plotted as a function of the 
ratio of blade chord to stall-zone width. As the stall-zone width in- 
creases (l.e., smaller values of c/z), the propagation rate increases, 
as Stenning predicted. Extrapolating the experimental curve to c/z = 0 
corresponds to the actuator disk Marble and Sears studied (one of infi- 
nitely small chord) . The values of stall-propagation rate obtained from 
the Sears theories are shown at c/z = 0 for comparison. The stall- 
propagation rates obtained from the Sears channel theory (eq^. (4)) with a 
phase shift 6 of 5° and P of 60° and from the airfoil theory (eq. (2)) 
with a phase shift A of 80° and p of 60° are in good agreement with 
the extrapolated data. A x value of 0.824 is required to make Stenning' s 
predicted propagation rate (eq. (5)) agree at c/z = 0 for p = 60° and 




SUMMARY OF RESULTS 

The rotating- stall characteristics of a 0.9 hub-tip ratio rotor were 
investigated, with the following results and observations : 

1. As the flow coefficient was reduced, total-span stall zones de- 
veloped in the following order: two zones, three zones, and one zone. 
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The one-stall-zone pattern, which produced the most severe pressure, tem- 
peratvtre, and flow fluctuations, spread over approximately one-third of 
the annulus . 

2. The rotor performance in the flow area between the stall zones 
was in good agreement with that obtained under steady-state conditions 
at the same flow. 

3. Static-pressure fluctuations downstream of the rotor were about 
60 percent as large as and approximately 180° out of phase with the 
pressure fluctuations upstream of the rotor. 

4 . The relative stall-propagation rates increased with increasing 
stall-zone size. 

5. The measured stall-propagation rates can be made to agree with 
theoretical predictions by adjusting certain theoretically undetermined 
factors appearing in the theories. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, Jvme 13, 1955 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

C 2 constants in equation for heat-transfer characteristics of hot- 
wire anemometer (eq. (B12)) 

blade chord 

specific heat at constant pressure 

voltage drop across hot-wire anemometer^ v 

voltage drop across hot-wire anemometer with zero flow, v 

indicates functional relation with quantity in parentheses 

frequency with which disturbance passes hot-wire anemometer or 
pressure transducer, cps 

acceleration due to gravity, ft/sec^ 

absolute rotative speed of a propagating-stall zone, rps 
current through hot-wire anemometer, ainp 
mechanical equivalent of heat, ft-lb/Btu 
integers (appendix B) 

slope of curve of rotor total -pressure-rise coefficient against 
inlet angle 

slope of curve of blade lift coefficient against inlet angle 
compressor rotative speed, rps or rpm as indicated 
static pressure, Ib/sq ft 
Ib/sq ft 

ratio of local radius to rotor tip radius 
total temperature , 
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u 

rotor blade speedy ft/sec 

V 

absolute velocity^ ft/sec 

w 

relative velocity, ft/sec 

X 

time between appearance of a flow fluctuation on oscillogram 

1 

* 

from one of two anemometers and appearance of a flow fluctu- 
ation on other anemometer (see fig. 23), sec 

z 

width of stall zone 

a 

angular displacement between two hot-wire anemometers, deg 

3 

relative inlet air angle measured from axial direction, deg 

r 

relative outlet air angle measured from axial direction, deg 

A 

denotes finite difference of quantity following 

A 

phase-shift angle (refs. 7 and 8), deg 

6 

phase-shift angle (ref. 8), deg 

e 

ratio of compressor -inlet temperature to NACA standard sea- 
level temperature of 518.70 R 

X 

effective blade passage area reduction (ref. 10) 

X 

number of stall zones in annulus 

V 

angular displacement between two propagating-stall regions, deg 

P 

density, slugs/cu ft 


mass -flow rate indicated by average current through hot wire, 
lb/(sec)(sq ft), (eq. (B17)) 

ApV/^ 

amplitude of flow fluctuation divided by pV 

a 

cascade solidity 

<(> 

flow coefficient, Vz/^m 
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^ pressure coefficient^ ^/(pU^n^/s) 

2 resistance of hot-wire anemometer^ ohms 

Subscripts : 

e effective gas 

i used to associate value of a and possible corresponding 

values of X 

m mean radius (r = 0.95) 

s stall 

w hot-wire anemometer 

z axial 

1 , 2 , 

. . .8 axial stations (fig. 4) 

Presubscripts : 

0 , 1 , 

. . .k indicates value of n used in determining X (see appendix B) 
Superscripts : 

+ indicates clockwise rotation of propagating stall 

indicates counterclockwise rotation of propagating stall 
^ indicates differentiation 

(bar) time-averaged value 
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APPENDIX B 

METHOD OF DETERMINING NUMBER OF STALL ZONES AND AMPLITUDE 
OF FLOW FLUCTUATIONS DUE TO ROTATING STALL 
Number of Stall Zones 

In order to derive the expression necessary to determine the num- 
ber of stalled regions by use of two anemometers displaced angularly in 
the plane of instrumentation (fig. 23(a)), it will be convenient to 
assume that the number of stalls and their direction of rotation are 
known . 

A stalled region passes each anemometer station / times per sec- 
ond, and 


/ = hX (Bl) 

The output signals of the two anemometers displaced angularly a de- 
grees appear on the oscilloscope as shown in figure 23(b) . If the 
stalled regions are rotating clockwise as indicated in figure 23(a), a 
stall zone will pass each anemometer every l// second, and a given 
stall zone will first pass anemometer 1 and then pass anemometer 2 at 
a/360h seconds later. Inasmuch as the stalls are quite similar, one 
stall zone cannot be distinguished from another, so a/360h cannot, 
in general, be determined from a single oscillogram. The phase shift 
(i.e., the time x) can be read directly from the oscillogram. For the 
example shown, 

X = - - 

36Ch / (B2) 

Then the ratio ^ - 1 = y. 

Since 

360 360h 

then 

y = ^ - 1 (B3) 

If there are twice as many stalled regions. 


360h 


2 

/ 
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Then, in general, with the stalls rotating in the direction shown in 
figure 23(a), 


and 


X 


a n 
360h ■ / 


y = f - n (B4) 

where n is the number of stalls enclosed in angle a and consequent- 
ly an integer. The value of y then varies with a/v as shown in 
figure 23(c) . When equation (B4) is differentiated with respect to a, 

^ i = -LL 

da V 360 

is obtained, provided a/^ is not an integer, in which case the deriva- 
tive does not exist in the usual sense because y is discontinuous at 
these points. 



If the stalls are rotating in the direction opposite (counter- 
clockwise) to that indicated in figure 23(a), a stall will pass ane- 
mometer 2 before passing anemometer 1. For counterclockwise rotation, 
the oscillograms would appear as shown in figure 23(d) . In this case. 


and 


X 


2 a 
7 ” 360h 



(B6) 


If there are twice as many stalls, 

y “ 3 - 55: 

V 


and, in general. 


y = n + 1 - 5!: 

V 


(B7) 


A plot of y against a/v is shown in figure 23(e), where n is, as 
before (eq. (B4) ) , the number of stalls enclosed in angle a and con- 
sequently an integer. Differentiation of (B7) with respect to a gives 

da V 360 
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provided a/v is not an integer, in which case the derivative (BS) does 
not exit. 

In order to determine the number of stalls X in any given case, 
two anemometers are needed with provision for varying the angle a. If 
one of the anemometers is on a movable mount such that a may be varied^ 
continuously, X can be determined unambiguously from equation (B5) or 
(B8). 


When the anemometers are located in fixed positions such that a 
is changed hy steps ^ the number of stalls n included in each value of 
a must be considered. From equation (B4) and 

V = 

X 


there is obtained 


_ 

^ " 360 


n 


(B9) 


or 


.+ , V 360 

X+ = (y -H n) — 


for stall zones rotating clockwise. From equation (B7), 


y = n + 1 


gX 

360 


or 


X" = Cn + 1 - y) ^ (BIO) 

ULr 

for counterclockwise rotation of the stalled regions. 

For various values of a, y is determined from the oscillograms; 
then X may be computed for values of n = 0, 1, 2,3,.. .k by use 
of equations (B9) and (BIO) and tabulated as shown in the following 
table. Unless y is measured quite accurately, X may not be an in- 
teger; however, the nearest integer should be tabulated. 
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a 

ap 

y 

yi 

n 

0 

1 

2 

3 

k 

x+ 

1 — 1 
o 

X+ 

1 1 

1 — 1 
00 

X+ 

3 1 

x+ 

k 1 

x" 

1 — 1 
o 

1 — 1 
1 — 1 

1 — 1 
00 

3^'l 

1 — 1 


^2 

ap 

yz 

yi 

0 

1 

2 

3 

k 

0 

1 

2 

3 

k 

X*” 

0^2 

1^"2 



.’'"a 

•H 

o 


E>^^ 

X+ 

3 i 

X+ 

k i 

0^2 



3^2 


0 i 


e’^'i 

X" 

3 i 

X“ 

k i 


The correct value of the number of stalls X and their direction 
of rotation will appear for every value of anemometer spacing a and 
the corresponding value of y measured from the oscillogram. If the 
rotation is clockwise, the correct value of X will recur in the 
column \ ; if the rotation is counterclockwise^ the correct value of 
\ will recur in the column X". If more than one value of X recurs, 
additional values of a must be tried to determine X* 

Perhaps a more convenient way of determining X is to plot equa- 
tions (B9) and (BIO) . A schematic sketch of such plots of X+ and X“ 
for a particular a is shown in figure 24. Several values of a may 
be included, but only one is shown here for illustrative purposes. By 
entering the ordinate at the y value determined from oscillograms, 
the X values, read to the nearest integer on the abscissa, occur at 
the intersections of the y value and the lines for the integer n. 

As described before with the tabular method^ the correct value for X 
will recur for every value of anemometer spacing a. 


Flow Fluctuations 

In addition to determining the number of stall zones, the hot-wire 
anemometer is also used to provide information regarding the magnitude 
of flow fluctuations during stall. The oscillogram of the hot -wire - 
anemometer output is a measure of the voltage drop across the wire 
element of the anemometer probe necessary to maintain it at a fixed 
resistance and consequently at a constant temperature. The time- 
averaged value of the anemometer -probe current is proportional to the 
time-averaged value of voltage drop across the wire element, 




e s= i 2 


(Bll) 
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The voltage drop e is related to the flow rate pV past the 
wire element by its heat-transfer characteristics. Reference 15 in- 
dicates that the heat-transfer characteristics of the wire element can 
be correlated satisfactorily by the following formula (King's law); 

1^2 = [C^ + f(M)C2V^j (T^ - Tg) (B12) 


The factor f(M) is a function of Mach number, and both Cj_ and 83 

are determined by calibration. Solving (B12) for pV results in the 
following expression: 


pV = 




-1 


C2f(MKTv - Tg) “ C2f(M) 
By use of the relation e = i2, 

pV = 


2 

e^ 1 


Si(T^ - Tg)f(M)C2 " f(M)C2 


(B13) 


(B14) 


If the effect of Mach number of f(M) and possible variation in Tg 
during a fluctuation in pV are neglected. 


pV = f(e) 


(B15) 


Fluctuations in pV will cause variations in e about some time- 
averaged value e; then 


pV = f(e) = f(e + Ae) 


(B16) 


Inasmuch as primary interest is in the magnitude of the pV fluc- 
tuations compared with the average value, it will be convenient to ex- 
press the pV fluctuations in terms of the voltage fluctuation across 
the wire element of the anemometer probe as measured from the oscillo- 
grams . Defining 


pV = 


_2(T^ - Te)C2f(M) ' ?(mJ^ 
and expanding (BlS) in a Taylor series give 


= f(e) 


f(e) . f(e) + Ae f(l) * , iM 


41 


(B17) 




(BIS) 
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Using 


where 


f(e) = V - B, 


A = 2(T^ - T^)C2f(M) 


(A and B are assumed constant during fluctuation in pV)^ 


1 

O ^ 


f(e) [f(a3^ 

1 

...(el - df(e)]^ ^ sf 

■ - A a2 




f " ' (e) = 


24e 


/-X 24 

f --'(e) = 




When (B20) is substituted in (B18) , 


1 

2 


f(e) = f(^) + ^2 


/Ae'\^ (t^\ 


+H-; *\~)j 


&(®)] « 


(B19) 


(B20) 


‘ f * <f)‘ 

(B21) 


results. Rearranging (B21) and dividf jig toth sides by f(e) give 

,2 /a^\3 /a^\ 4 


ne) - [y, 

A^ffe) L \ 


A‘^f(e) 
2 

T 

2 


Ae 


<f) -(f). 

t f - ^ (f )!1 


(B22) 


A[f(e)] 
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From (B16) and (B17), 


The factor 


f(e) - f(e) _ pV - pV ^ ApV 
f(e) ^ ^ 


e2 e2 1 


A[f(e)] 


1 - /-2 \ " e 2 - AB 1 - M 

e2 


~ M 2 Al 


- b ) 


From (B19 ) , 


^ ^ CiKTy - Te) 

^ ~ e2 


(B23) 


(B24) 


but 


When pV = 0^ e = eQ^ so 


and 


C-i = 


^0 


1 2(T^ - Tj 


i - (?/ 




A[f(e)]5 l-(^f 

Substituting (B25) and (B23) in (B22) yields 


(B25) 


AdV ^ 


i-(¥ 


'4 f 4 4 


(f J - (f )' 


-(¥/ 


+ 


(B26) 


A plot of A^V/pV against Ae/e for various values of Bq/s is pre- 
sented in figure 25. The mass-flow rate pV is zero when 


Ae 


e er 


e 


e 
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The frequency with which the stall zones passed the hot wire was 
determined by means of Lissajous figures. Such figures can be formed 
by putting the hot-wire signal on the y-axis of the oscilloscope and 
an audio -frequency oscillator on the x-axis. The oscillator is 
adjusted until the trace on the oscilloscope appears to be a stationary 
near-circle. 
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(a) Discontinuous. 



Flow coefficient, V 
(b) Continuous. 

Figure 1. - Compressor performance characteristics. 
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Figure 2. - Vector velocity diagram for 0.9 hub-tip 
ratio rotor at r = 0.95. 
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Figure 3. - Schematic diagram of con 5 )ressor installation. 
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Figure 4. - Axial measuring stations along rotor system. 
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Figure 5. - Installation of pressure pickup in compressor casing. 
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Figure 6. - Static performance characteristics of rotor 
test system. 
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Figure 7 . - Rotor 
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at r = 0.95. 
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Figure 8. - Flow fluctuations along rotor unit axis. 
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Time ^ 

(a) Two stall zones; stall frequency, 175 cps; flow coefficient, 0.75. 

Figure 9. - Hot-wire-anemometer oscillograms at stations 2, 4, 7^ and 8 for 
0.9 hub- tip ratio rotor. Typical axial sur\^ey at r = 0.95 for rotor speed 
of 8000 rpm. 
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Station Amplitude, 

AoV/pV 

4 0.13 


2 .03 


4 ,10 


7 .36 


4 .14 


C- 38984 


Time — 

(b) Three stall zones; stall frequency, 2J6 cps; flow coefficient, 0.68. 

Figure 9. - Continued. Hot-wire- anemometer oscillograms at stations 2, 4, 7, 
and 8 for 0.9 hub-tip ratio rotor. Typical axial sirrvey at r = 0.95 for 
rotor speed of 8000 rpm. 


8 .25 






36 


MCA TN 3518 




4 


7r 


4 


8 

C- 38985 


1.17 


1.12 


‘ 1.19 


1.18 


Time — 

(c) One stall zone; stall frequency^ 60 cps; flow coefficient, 0.60. 

Figure 9. - Concluded. Hot-wire anemometer oscillograms at stations 2, 4, 7, 
and 8 for 0.9 hub- tip ratio rotor. Typical axial survey at r = 0.95 for 
rotor speed of 8000 rpm. 
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Flow-fluctuation amplitude, ZipV/pV 


Figure 10. - Flow fluctuations along blade span at station 5. 
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ratio, r 
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Figure 11. - Hot-vire- anemometer oscillograms at stations 4 and 5 for 0.9 
hub- tip ratio rotor. Typical radial survey at rotor speed of 8000 rpm. 
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(Td) Three stall zones; stall frequency, 276 cps; flow coefficient, 0.68. 

Figure 11. - Continued. Hot-wire- anemometer oscillograms at stations 4 and 5 
for 0.9 hub-tip ratio rotor. Typical radial survey at rotor speed of 
8000 rpm. 
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Station Amplitude , Eadius 

Apv/pV ratio, r 


4 1.09 0.95 


5 1.01 .98 


4 1.06 .95 


5 1.07 .95 




Time — ^ 

(c) One stall zone; stall frequency, 60 cps; flow coefficient, 0.60. 

Figure 11. - Concluded. Hot-wire-anemometer oscillograms at stations 4 and 5 
for 0.9 hub-tip ratio rotor. Typical radial survey at rotor speed of 
8000 rpm . 
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(a) Two stall zones; stall frequency^ 1J5 cps; flow coefficient^ 0.705. 

Figure 12. - Oscillograms from high- frequency-response pressure pickups 
for 0.9 huh-tip ratio rotor at r = 1.0 and rotor speed of 8000 i*pm. 
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Time — ► 

("b) Three stall zones; stall frequency^ 276 cps; flow coefficient, 0,68. 

'Fl^are 12. - Continued. Oscillograms from hi^-frequency-response pressure 
pickups for 0.9 hub-tip ratio rotor at r = 1.0 and rotor speed of 8000 rpm. 
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Time — ^ 

(c) One stall zone; stall frequency^ 60 cps; flow coefficient, 0.63. 

Figiare 12. - Concluded. Oscillograms from high-frequency-response pressure 
pickups for 0.9 huh-tip ratio rotor at r = 1.0 and rotor speed of 8000 rpm. 
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Station 3 


C-38992 


Time — ^ 

Figure 13, - Phase relation hetween pressure fluctuations at stations 
7 and 3 for 0.9 hub- tip ratio rotor at r = 1.0 and rotor speed of 
8000 rpm. One stall zone; stall frequency, 60 cps; flow coefficient, 
0.63. 
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(a) Two stall zones; flow coefficient^ 0.705. 


Figure 14. - Sketch from oscillogram of static-pressure variation 
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Time 

(b) Three stall zones; flow coefficient, 0.68. 


Figure 14. - Continued. Sketch from oscillogram of static- 
pressure variation. 
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(c) One stall zone; flow coefficient, 0.63. 

Figure 14. Concluded. Sketch from oscillogram of static-pressure variation. 
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Pressure coefficient across rotor, ^ 3.7 
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Figure 15. - Variation of pressure coefficient across 
rotor with time. One stall zonej flov coefficient, 
0.63. 
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Flov coefficient, 4> 


Figure 16. - Static performance characteristics between 
stations 1 to 3 and 3 to 7 . 
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Flow coefficient, (p 


Figure 17. - Average temperature rise measured across rotor. 
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Time 

Figure 18. - Variation of temperature rise with time for single-zone stall 
pattern with flow coefficient of 0.63. 
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. - Stall propagation in cascade. 


Figure 19 






Figure 20. - Comparison of Sears airfoil theory (ref. 7) and experimental 
absolute stall-propagation rates. 
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Relative inlet-air angle at stall, p, deg 


Figure 21. - Comparison of Sears channel theory (ref. 8) and Marble theory 
(ref. 9) with experimental absolute stall -propagation rates. 
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0 «05 • 10 *15 *20 «25 

Ratio of "blade chord to stall-zone width, c/z 

Figure 22. - Effect of stall-zone size on relative 
stall -propagation rate. 
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Anemometers 

1 2 0 <a< n 



(a) Sketch showing anemometer location and propagating 
stall in compressor annulus. 



(h) Oscillogram from cathode-ray oscilloscope, clockwise rotation of 

stall. 



(c) Variation of y with a/^ , clockwise rotation of stall. 

Figure 23. - Relations necessary for determining number of 
propagating stalls. 
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(d) Oscillogram from cathode-ray oscilloscope, counter- 
clockwise rotation of stall. 



a/ V 

(e) Variation of y with cl/v, counterclockwise rotation 

of stall. 

Figure 23. - Concluded. Relations- necessary for determining 
number of propagating stalls. 
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(Td) CounterclockwiBe rotation of stall zone. 

Figure 24. - Schematic representation of plots useful in determining 
number of stall zones. 


Plow fluctuation, ApV/^ 


NACA TW 3518 


59 



-.5 -.4 -.3 -.2 -.1 0 .1 .2 .3 

Anemometer voltage fluctuation, Ae/e 

Figure 25. - Variation of weight-flow fluctuations with anemometer 
voltage fluctuation. 
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